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Abstract 

A quantum control landscape is defined as the physical objective as a function of 
the control variables. In this paper the control landscapes for two-level open quantum 
systems, whose evolution is described by general completely positive trace preserving 
maps (i.e., Kraus maps), are investigated in details. The objective function, which is 
the expectation value of a target system operator, is defined on the Stiefel manifold 
representing the space of Kraus maps. Three practically important properties of the 
objective function are found: (a) the absence of local maxima or minima (i.e., false 
traps); (b) the existence of multi-dimensional sub-manifolds of optimal solutions corre- 
sponding to the global maximum and minimum; and (c) the connectivity of each level 
set. All of the critical values and their associated critical sub-manifolds are explicitly 
found for any initial system state. Away from the absolute extrema there are no local 
maxima or minima, and only saddles may exist, whose number and the explicit struc- 
ture of the corresponding critical sub-manifolds are determined by the initial system 
state. There are no saddles for pure initial states, one saddle for a completely mixed 
initial state, and two saddles for partially mixed initial states. In general, the landscape 
analysis of critical points and optimal manifolds is relevant to explain the relative ease 
of obtaining good optimal control outcomes in the laboratory, even in the presence of 
the environment. 

1 Introduction 

A common goal in quantum control is to maximize the expectation value of a given target 
operator by applying a suitable external action to the system. Such an external action 
often can be realized by a tailored coherent control field steering the system from the initial 
state to a target state, which maximizes the expectation value of the target operator [U El 
El m El El [71 El |9]. Tailored coherent fields allow for controlling Hamiltonian aspects (i.e., 
unitary dynamics) of the system evolution. Another form of action on the system could 
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be realized by tailoring the environment (e.g., incoherent radiation, or a gas of electrons, 
atoms, or molecules) to induce control through non-unitary system dynamics [10]. In this 
approach the control is the suitably optimized, generally non-equilibrium and time dependent 
distribution function of the environment; the optimization of the environment would itself 
be attained by application of a proper external action. Combining such incoherent control 
by the environment (ICE) with a tailored coherent control field provides a general tool 
for manipulating both the Hamiltonian and dissipative aspects of the system dynamics. A 
similar approach to incoherent control was also suggested in [11] where, in difference with [TO], 
finite-level ancilla systems are used as the control environment. The initial state of the field 
and the interaction Hamiltonian as the parameters for controlling non-unitary dynamics 
was also suggested in [12]. Non-unitary controlled quantum dynamics can also be realized 
by using as an external action suitably optimized quantum measurements which drive the 
system towards the desired control goal [TSj [HI [151 [IS [H] • General mathematical definitions 
for the controlled Markov dynamics of quantum-mechanical systems are formulated in |18j . 

In this paper we consider the most general physically allowed transformations of states 
of quantum open systems, which are represented by completely positive trace preserving 
maps (i.e., Kraus maps) [TH [20l [211 122]. A typical control problem in this framework is 
to find, for a given initial state of the system, a Kraus map which transforms the initial 
state into the state maximizing the expected value (O) of a target operator of the system. 
Practical means to find such optimal Kraus maps in the laboratory could employ various 
procedures such as adaptive learning algorithms [S, [23], which are capable of finding an 
optimal solution without detailed knowledge of the dynamics of the system. Kraus maps 
can be represented by matrices satisfying an orthogonality constraint (see Sec. II), which can 
be naturally parameterized by points in a Stiefel manifold [21], and then various algorithms 
may be applied to perform optimization over the Stiefel manifold (e.g., steepest descent, 
Newton methods, etc. adapted for optimization over Stiefel manifolds) [251 [26] . 

The quantum control landscape is defined as the objective expectation value (6) as a 
function of the control variables. The efficiency of various search algorithms (i.e., employed 
either directly in the laboratory or in numerical simulations) for finding the minimum or 
maximum of a specific objective function can depend on the existence and nature of the 
landscape critical points. For example, the presence of many local minima or maxima (i.e., 
false traps) could result in either permanent trapping of the search or possibly dwelling for a 
long time in some of them (i.e., assuming that the algorithm has the capability of extricating 
the search from a trap) thus lowering the search efficiency. In such cases stopping of an 
algorithm at some solution does not guarantee that this solution is a global optimum, as 
the algorithm can end the search at a local maximum of the objective function. A priori 
information about absence of local maxima could be very helpful in such cases to guarantee 
that the search will be stopped only at a global optimum solution. This situation makes 
important the investigation of the critical points of the control landscapes. Also, in the 
laboratory, evidence shows that it is relatively easy to find optimal solutions, even in the 
presence of an environment. Explanation of this fact similarly can be related with the 
structure of the control landscapes for open quantum systems. 

The critical points of the landscapes for closed quantum systems controlled by unitary 
evolution were investigated in [23 [23 [23 [30l [31] , where it was found that there are no sub- 
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optimal local maxima or minima and only saddles may exist in addition to the global maxima 
and minima. In particular, it was found that for a two-level system prepared initially in a 
pure state the landscape of the unitary control does not have critical points except for global 
minima and maxima. 

The capabilities of unitary control to maximize or minimize the expectation value of the 
target operator in the case of mixed initial states are limited, since unitary transformations 
can only connect states (i.e., density matrices) with the same spectrum. In going beyond 
the latter limitations, the dynamics may be extended to encompass non-unitary evolution by 
directing the controls to include the set of Kraus maps (i.e., dual manipulation of the system 
and the environment). Quantum systems which admit arbitrary Kraus map dynamics are 
completely controllable, since for any pair of states there exists a Kraus map which transforms 
one into the another [32] . 

In this paper the analysis of the landscape critical points is performed for two-level 
quantum systems controlled by Kraus maps. It is found that the objective function does 
not have sub-optimal local maxima or minima and only saddles may exist. The number of 
different saddle values and the structure of the corresponding critical sub-manifolds depend 
on the system initial state. For pure initial states the landscape has no saddles; for a 
completely mixed initial state the landscape has one saddle value; for other (i.e., partially 
mixed) initial states the landscape has two saddle values. For each case we explicitly find all 
critical sub-manifolds and critical values of the objective as functions of the Stokes vector of 
the initial density matrix. An investigation of the landscapes for multi-level open quantum 
systems with a different method may also be performed [36|. The absence of local minima or 
maxima holds also in the general case although an explicit description of the critical manifolds 
is difficult to provide for multi-level systems. The absence of false traps practically implies 
the relative ease of obtaining good optimal solutions using various search algorithms in the 
laboratory, even in the presence of an environment. 

It should be noted that the property of there being no false traps relies on the assumption 
of the full controllability of the system, i.e., assuming that an arbitrary Kraus map can be 
realized. Restrictions on the set of available Kraus maps can result in the appearance of 
false traps thus creating difficulties in the search for optimal solutions. Thus, it is important 
to consider possible methods for engineering arbitrary Kraus type evolution of a controlled 
system. One method is to put the system in contact with an ancilla and implement, on 
the coupled system, specific unitary evolution whose form is determined by the structure 
of the desired Kraus map [37] (see also Sec. II). Lloyd and Viola proposed another method 
of engineering arbitrary Kraus maps, based on the combination of coherent control and 
measurements [38j. They show that the ability to perform a simple single measurement on 
the system together with the ability to apply coherent control to feedback the measurement 
results allows for enacting arbitrary Kraus map evolution at a finite time. 

A level set of the objective function is defined as the set of controls which produce the 
same outcome value for (0). We investigate connectivity of the level sets of the objective 
functions for open quantum systems and show that each level set is connected, including the 
one which corresponds to the global maximum/minimum of the objective function. Connec- 
tivity of a level set implies that any two solutions from the same level set can be continuously 
mapped one into another via a pathway entirely passing through this level set. The proof of 
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Figure 1: This figure schematically illustrates the landscape J as a function of two controls xi 
and X2- The figure shows the two main properties of quantum-mechanical control landscapes 
for open quantum systems: (a) absence of false traps and (b) connectivity of the sub-manifold 
of global maximum solutions (a one dimensional curve at the top of the landscape in this 
example) . 

the connectivity of the level sets is based on a generalization of Morse theory. Experimental 
observations of level sets for quantum control landscapes can be practically performed, as it 
was recently demonstrated for control of nonresonant two-photon excitations [39j . 

In summary, the main properties of control landscapes for open quantum systems are: 
(a) the absence of false traps; (b) the existence of multi-dimensional sub-manifolds of global 
optimum solutions, and (c) the connectivity of each level set. The proof of the properties (a)- 
(c) is provided in the next sections for the two-level case. Figure [l] illustrates the properties 
(a), (b), and connectivity of the manifold of global maximum solutions; the figure does 
not serve to illustrate other properties such as connectivity of each level set. It is evident 
that the function drawn on figure [T] does not have local minima or maxima and the set of 
solutions for the global maximum is a connected sub- manifold (a curve in this case). A 
simple illustration is chosen for the figure since an exact objective function for an A^-level 
quantum system depends on D = 2A^^ — 2N'^ real variables (such that D = 24 for N = 2) 
and therefore can not be drawn. 

The present analysis is performed in the kinematic picture which uses Kraus maps to 
represent evolution of quantum open systems. An important future task is to investigate 
the structure of the control landscape in the dynamical picture, which can be based on 
the use of various dynamical master equations to describe the dynamics of quantum open 
systems [221 SHI IHl 1121 SSI SI] ■ Such analysis may reveal landscape properties for quantum 
open systems under (possibly, restricted) control through manipulation by a specific type of 
the environment (e.g., incoherent radiation). 

In addition to optimizing expected value of a target operator, a large class of quantum 
control problems includes generation of a predefined unitary (e.g., phase or Hadamard) [21j or 
a non-unitary [33] quantum gate (i.e., a quantum operation). This class of control problems 
is important for quantum computation and in this regard a numerical analysis of the problem 
of optimal controlled generation of unitary quantum gates for two-level quantum systems 
interacting with an environment is available [3ll 135] . 
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Although the assumption of complete positivity of the dynamics of open quantum sys- 
tems used in the present analysis is a generally accepted requirement, some works consider 
dynamics of a more general form [ISl HH] • Such more general evolutions may result in dif- 
ferent controllability and landscape properties. For example, for a two-level open quantum 
system positive and completely positive dynamics may have different accessibility proper- 
ties [47j. In this regard it would be interesting to investigate if such different types of the 
dynamics have distinct essential landscape properties. 

In Sec. [2] the optimal control problem for a general A^-level open quantum system is 
formulated. Section |3] reduces the consideration to the case of a two-level system. In Sec.|4]a 
complete description is given of all critical points of the control landscape. The connectivity 
of the level sets is investigated in Sec. [5j 

2 Formulation for an A^- level system 

Let AiN be the linear space oi N x N complex matrices. The density matrix p of an A^-level 
quantum system is a positive component in A^at, p > 0, with unit trace, Trp = 1 (Hermicity 
of p follows from its positivity). Physically allowed evolution transformations of density 
matrices are given by completely positive trace preserving maps (i.e., Kraus maps) in Ai^. 
A linear Kraus map ^ : A4n —y A4n satisfies the following conditions [191 '■ 

• Complete positivity. Let I„ be the identity matrix in Ain- Complete positivity means 
that for any integer n G N the map $ ® In acting in the space A^at ® Ain is positive. 

• Trace preserving: Vp G A^at, Tr$(p) = Trp. 

Any Kraus map $ can be decomposed (non- uniquely) in the Kraus form [9] : 

M 



where the Kraus operators Ki satisfy the relation Yli=i ^i^i ~ ^n- For an iV-level quantum 
system it is sufficient to consider at most M = N"^ Kraus operators |48j . 

Let Til = be the Hilbert space of the system under control. An arbitrary Kraus map 
of the form ([T]) can be realized by coupling the system to an ancilla system characterized 
by the Hilbert space 7^2 = C^, and generating a unitary evolution operator U acting in 
the Hilbert space of the total system H = Hi ® 0.2 as follows [37j. Choose in TC2 a unit 
vector |0) and an orthonormal basis |ej), i = 1, . . . , M. For any \ip) G Hi let U^ip) ® |0)) = 
^^iKi\ip) ® \ei). Such an operator can be extended to a unitary operator in H and for 
any p one has $(p) = Ty-h^ {U{p ® |0)(0|)[/"'"}. Therefore the ability to dynamically create, 
for example via coherent control, an arbitrary unitary evolution of the system and ancilla 
allows for generating arbitrary Kraus maps of the controlled system. 

Let po be the initial system density matrix. A typical optimization goal in quantum 
control is to maximize the expectation value J = (6) of a target Hermitian operator O over 
an admissible set of dynamical transformations of the system density matrices. For coherent 
unitary control this expectation value becomes 




(1) 



1=1 



J[U] = Ti[UpoU^Q] 
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where U = U{t,to) is a unitary matrix, = = In, which describes the evolution 
of the system during the control period from the initial time to until some final time t and 
implicitly incorporates the action of the coherent control field on the system. 

In the present paper we consider general non-unitary controlled dynamics such that the 
controls are Kraus maps, for which the parametrization by Kraus operators is used. The 
corresponding objective function specifying the control landscape has the form 



where the Kraus operators {Ki} = {Ki{t, to)} describe evolution of the open quantum system 
from an initial time to until some final time t. The control goal is to maximize the objective 
function over the set of all Kraus operators Ki, . . . , Km satisfying Xlfii ^l^i — ^n, thereby 
forming a constrained optimization problem. 

Definition 1 Let ¥ be a field of real or complex numbers, i.e., ¥ = M. or ¥ = C A Stiefel 
manifold over ¥, denoted 14 (F"), is the set of all orthonormal k-frames in F" (i.e., the set 
of ordered k-tuples of orthonormal vectors in ¥"-). The case F = M (respectively, ¥ = C) 
corresponds to a real (complex) Stiefel manifold. 

Let K be the Nx{NM) matrix defined as i^' = [K^ . . . Kj^), where is the transpose of 
matrix Ki and M is the number of Kraus operators. Consider vectors Xi, . . . , € 
with components {Xi)j = Kij, i.e., vector Xi is the z-th row of the matrix K. The constraint 
SiHi ^i^i = Iat in terms of the vectors Xi, . . . , X^ takes the form {Xi, Xj) = 6ij, where Sij is 
the Kronecker delta symbol. This constraint defines the complex Stiefel manifold V/v(C^^^). 
Therefore optimization of the objective function J[Ki, . . . ,Km] defined by Eq. ^ can be 
formulated as optimization over the complex Stiefel manifold Viv(C^*^). 

3 Two- level system 

In the following we consider the case of a two-level system in detail. Any density matrix of 
a two-level system can be represented as 



where a = (cxi, (T2, cts) = {ax, cry, az) is the vector of Pauli matrices and w G M'^ is the Stokes 
vector, ||w|| < 1. Thus, the set of density matrices can be identified with the unit ball in 
M^, which is known as the Bloch sphere. 

Any Kraus map $ on M.2 can be represented using at most four Kraus operators 



M 




(2) 



1=1 



p=-[l + (w,a)] 




/ = 1,2,3,4 



as $(p) = Ylii=i^iP^i 1 where the Kraus operators satisfy the constraint 
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(3) 



1=1 
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Let po be the initial system density matrix with Stokes vector w = {a,f3,'j), where 
||w|p = + + 7^ < 1, and let 6 be a Hermitian target operator. The objective 
functional for optimizing the expectation value of 9 has the form J[Ki, K2, K^, K^, po, ©] = 
Ylt=i Tr[-f^/Po-^/0] • The control goal is to find all quadruples of Kraus operators {Ki, K2, K^, K4) 
which maximize (or minimize, depending on the control goal) the objective functional J. 
The goal of the landscape analysis is to characterize all critical points of J[Ki, K2, K^, ^^4], 
including local extrema, if they exist. 

The analysis for an arbitrary 2x2 Hermitian matrix G can be reduced to the case 

which we will consider in the sequel. This point follows, as an arbitrary Hermitian operator 
O G A^2 has two eigenvalues Ai and A2 and can be represented in the basis of its eigenvectors 
as 

where Ai > A2. One has 9 = (Ai — A2)6o + A2I2 and 

4 

1=1 

4 4 

= (Ai - A2) J2 Tr[i^zPo^/eo] + Tr[^/Po^/] 
1=1 1=1 

= {Xi-X2)J[KuK2,Ks,K^;po,eo] + \2 

Therefore, the objective function for a general observable operator O depends linearly on the 
objective function defined for Qq. We denote J[Ki, K2, K3, K^, w] := J[Ki, K2, K^, po? ©o]- 
In the trivial case G = I2 the landscape is completely flat and no further analysis is needed. 



4 The critical points of the objective function land- 
scape 

The Kraus operators for a two-level system can be parameterized by a pair of vectors X, F G 
= (j^4 q£ form X = ui ® vi and Y = U2 ® V2, where ui = (xn, 0:21, xai, 0:41), 
vi = (a;i2,a;22,a;32,a;42), U2 = (xis, X23, 3:33, 2:43), and V2 = (xi4, X24, 3:34, 3:44) • The objective 
function in terms of these vectors has the form 

J[ui,U2,vi,V2;w] = ^ (l + 7)||nif + (l-7)||M2f + 2Re[2;o(ui,M2)] (4) 

where zq = a — i(3, {■,■) and || • || denote the standard inner product and the norm in 
(here the numbers a, f3, 7 are the components of the Stokes vector w = (a, f3, 7) of the initial 
density matrix po, see Sec. [3|. The constraint (Isj) in terms of the vectors X and Y has the 
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form ||X|| = = 1, {X,Y) = and determines the Stiefel manifold A4 = V2(C®). The 
matrix constraint ^ in terms of the vectors Ui and Vi has the form 



^3{ui,U2,Vi,V2) 



Ikilr + ||filr - 1 = (5) 

\\u2f + \\v2f-l=0 (6) 

{U1,U2) + {vi,V2) = (7) 



If zq 7^ 0, then the objective function is diagonahzed by introducing new coordinates 
(ui, U2, vi, V2) in C"*^^ according to the formulas 

ui = nui — iyu2, U2 = T-^uui + T-^fiU2 (8) 

Fo| l^^ol 

z* z* _ 

Vi = HVi - UV2, V2 = Y^t^Vi + 7-^MV2 (9) 

I I ^0 1 

where /i = \zq\/ a/2||w|| (||w|| — 7) and u = |zo|/\/2||w|| (||w|| +7). The objective function 
in these coordinates has the form 

J[x; w] = A+||Mif + A_||n2ir (10) 

where x = {ui,U2,Vi,V2) ^ -M. and A-t = (1 ± ||w||)/2. If zq = and 7 > (resp., 7 < 0), 
then the objective function ^ has the form (10) with Ui = Ui,Vi = Vi for i = 1,2 (resp., 



Ml = U2,U2 = Ui,Vi = V2,V2 = Vi). The constraints (|5])-([7]) in the new coordinates have the 
same form ^i{ui,U2, Vi, V2) = for z = 1, 2, 3. 

Theorem 1 Let w = (a,/?, 7) & M.^ be a real vector such that ||w|| < 1 and let X± = 
(1 ± ||w||)/2. For any such w, the global maximum and minimum values of the objective 
function J[ui, U2, Vi, V2] w] = A+||?iip + A_ \\u2\\'^ are 

min J[ui,U2,Vi,V2;w] = 

max J[ui,U2,vi,V2]Vf] = 1. 

The critical sub-manifolds and other critical values of J in are the following: 

Case 1. w = (the completely mixed initial state). The global minimum sub-manifold 
is A^min*^^ = {x ^ Ai \ Ui = U2 = 0} ■ The global maximum sub-manifold is A^max^"* = {x G 
I = ?)2 = 0}. The objective function has one saddle value J = 1/2 with the correspond- 
ing critical sub-manifold A^sadd^i = {x & M.\u2 = zui, Vi = —z*V2, 2; G C} \J{x G Ai | -Si = 
^2 = 0}. The Hessian of J at any point at M.f^^f]l has = 6 positive, i/_ = 6 negative, 
and vq = 16 zero eigenvalues. 

Case 2. < ||w|| < 1 (a partially mixed initial state). The global minimum sub-manifold 
is Aimin = {x & Ai\ui = U2 = 0}. The global maximum sub-manifold is Al^ax = {2; G 
M. \ vi = V2 = 0} . The objective function has two saddle values: 

, , 1 ± ||w|| , , , 

J±(w) = ^ = A±. (11) 
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The corresponding critical sub-manifolds are Ai"!^ = {x E Ai \ ui = V2 = 0} and Ai^ = {x G 
Ai \ U2 = vi = 0} . The Hessian of J at any point at Ai"^ (resp., AV^) has z/+ = 8 positive, 
z/_ = 6 negative (resp., = 6 positive, z/_ = 8 negative), and uq = 14 zero eigenvalues. 

Case 3. ||w|| = 1 (a pure initial state). The global minimum sub-manifold is A^mm ~ 
{x E Ai\ui = 0}. The global maximum sub-manifold is A^^ax = {a; G | = 0}. The 
objective function has no saddles. 

Proof. The objective function has the form J = pu, where pn is the diagonal matrix 
element of the density matrix. Therefore < J < 1 and the value J = (resp., J = 1) 
corresponds to the global minimum (resp., maximum). 



The constraints can be included in the objective function (10) by adding the term 
V, f]] = 77i$i + ri2^2 + 2Re [r^g^s], where the two real and one complex Lagrange multi- 
pliers r]i,ri2, and rj^ correspond to the two real and one complex valued constraints $1, $2, 
and $3, respectively. Critical points of the function J on the manifold Ai are given 
by the solutions of the following Euler-Lagrange equations for the functional J[u,v,X\ = 
J[u, v] + V, rj]: 



= Vii.J^ = {X+ + + V3U2 

= V^.J^ = r/3*Mi + (A_ + 7/2)^2 

= Vv* J ^ = r]ivi + r]3V2 

= Vs. J ^ = r]lvi + r]2V2 



(12) 
(13) 
(14) 
(15) 



where tti, -02, '^'i, "^2 satisfy the constraints ([5 )-([7| ). The proof of the theorem is based on the 
straightforward solution of the system ( 12 )-( 15 ). The case 2 will be considered first, followed 
by the cases 1 and 3. 

Case 2. < ||w|| < 1. Consider in M the open subset Oi = {a; G | ?)i 7^ 0, C2 7^ 0}. 
Let us prove that the set of all critical points of J in Oi is the set of all points of A\ such 
that Ml = M2 = 0. 

Suppose that there are critical points in Oi such that £ti 7^ or -112 7^ 0. For such points 
the following identity holds 



1^3 I 



(A+ + r/i)(A„ +7/2). 



In Oi, 7^ and therefore |r73| = 771772. This equality together with (16) gives 



V2 



-A, 



A+ 



(16) 



(17) 



Suppose that 773 7^ 0. Then, using (|12|) and (|14|), the constraint $3 gives 

0. 



(A+ + 7/1) II 771 1 



?7iii^'ir 



Constraint $1 gives \\Vi\\^ = 1 — Hwilp, and therefore rji = — A+Httip. Similarly we find 
r]2 = —A_ 11-0211^. Substituting these expressions for 7/1 and 772 into the (12) and (13) we find 



a/aIa^ 
v/AlT 



Ul\\\\U2W 
+ ||'Sl|P||772| 



A+(l 
A_(l 



\Ul\ 
l«2| 



\Ul\ 



^A_ 



A+||7ii||^ + ^/X^\\u2 

A+IImiIP + x/\^\\u2 
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This system of equations implies A_ = A+ -v^ w = which is in contradiction with the 



assumption ||w|| > for the present case. If 773 = 0, then it follows from (14), (15) that 



Vi = V2 = 0. In this case equations (12) and (13) have only the solution ui= 112 = 0. 

Points in Oi with Ui = U2 = Q form the global minimum manifold A^^in ~ ^2(C^), 
which is a Stiefel manifold and hence is connected. In some small neighborhood of zero 
we can choose ui and U2 as normal coordinates. So A^^in is non degenerate. Similar 
treatment of the region O2 = {x & M.\ui ^ U2 ^ ^} gives the global maximum manifold 
MZ.. = {xeM\v^=v2 = Q}. 

Now consider the region = {a; G | ^2 7^ 0, f 1 7^ 0}. In this region the objective 
function J has the form 

J[ui,U2.,vi,V2] = A_ + A+||mi||^ - A_||'U2||^. (19) 

Using the analysis for the region Oi, we conclude that the objective function has no critical 
points such that 7)2 7^ in Os. Therefore all critical points in are in the sub- manifold 
A/" = {x e At I ^2 = 0} C A^. The restriction of J to A^ has the form 

J[Ui,U2,Vi,V2]\m = A_ + A+||Mif . 

Note that A" is a subset of all sets of vectors {ui,U2,vi) satisfying the constraints 

||M2r = 1, ll""!!!^ + ll'^llP = I7 (Mi,M2)=0. 

It is clear from this representation of M that 'WJ\j\f = if and only if ui = 0. This gives the 
critical sub-manifold Al^ = {x & A4 \ Ui = V2 = 0}. The objective function has the value 
J |xw= A_ on this manifold. 

To show that this is a saddle manifold, and not a local maximum or minimum, we 
calculate the Morse indices of the objective function on Ai"^ and show that both positive and 
negative Morse indices are different from zero (the Morse indices are the numbers of positive, 
negative and zero eigenvalues of the Hessian of J and positive and negative Morse indices 
determine the number of local coordinates along which the function increases or decreases, 
respectively). With regard to this goal, consider the manifold /C := {x G C^^ | ^i{u,v) = 
0, $2(m, w) = 0}. Let X e At. Below we introduce some coordinates in a neighborhood of x 
on /C. 

For any z G such that z ^ we define the unit vector g{z) = z/\\z\\ G C^. Let 
(fi, i = 1, ... ,7 he some coordinate system on S''' (embedded in as a unit sphere with the 
origin at zero) in some neighborhood \4 of g{u2{x)) and ipi, i = 1, ... ,7 he some coordinate 
system on S*^ in some neighborhood Vy of g{vi{x)). We will use the following functions 
defined in some neighborhood of x on /C (2; G /C): 

(pi{z) = (piOgou2{z), i = l,...,7 
xjji{z) = ipiOgovi{z), i = l,...,7. 

Let TzS'^ he the maximal complex subspace of the tangent space of S*^. For each z G \4 
let xi, . . . , xe be coordinates on T^S"^ and for each 2; G K Z/i, • • • , 1/6 be coordinates on TzS"^ . 
Let xi, . . . , xe and yi, . . . ,yQ he functions on /C defined as follows. 
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Let z = {ui,U2,vi,V2) G /C be in a small enough neighborhood of x. By definition Pr^ 
is the projection from C*^ to Tg(^u2)S'^ ^-iid Pr^ is the projection from to Tg(^y^^S'^. By 
definition 



Xi 



Xi o Pr„ o ui, 
Hi o Pr^, o V2, 



,6, 
,6. 



Now let Pr'^ and Pr^ be the complex-valued functions defined on by the formulas 



Pr'n(/) 
Pr'.(/) 



{9{u2),f), 

{9(yi),f), 



/ G 



By definition 



p:=Pr'„oui, q ■=Pt'^ov2. 



Thus, the functions ipi, ipi, Xk, yi,p, q, where i,j = l,...,7 and k,l = 1, . . . ,6, are coordinates 
on /C in some neighborhood of the point x. Locally the manifold is a sub-manifold of /C 
defined by the constraint $3 = 0. In our coordinates this constraint has a form 



P 



i=l 



i=l 



Therefore (fi, ipi, Xk, yi,p, where i,j = 1, . . . , 7 and k,l = 1, . . . ,6 are the coordinates on Ai 
in some neighborhood of x. The second differential of J at the point x in this coordinates 
has the form 



dV = A_ 



6 

E 

i=l 



dxf - A_ 



6 

J2dyl + {K-)^-)\dp\'. 



Since A+ — A_ = ||w|| > for the present case, the Morse indices of this point are z/+ = 
8, z/_ = 6 (note that p is a complex coordinate). 

Similar treatment of the region O4, = {x & Ai \ ui ^ 0, V2 ^ 0} shows the existence of 
the critical sub-manifold = {x & Ai \ U2 = vi = 0}. This sub-manifold corresponds to 



the critical value J 



A+ and its Morse indices are = 6, 



this concludes the proof for the case < ||w|| < 1. 
Case 1. w=0. Consider in A4 the open subset Oi. 



8. Since [j d = M, 

i=0 



Let ?73 = 0. Then in the region Oi Eqs. (14) and (15) imply that rjiVi = 772^2 = ^ rji = 
ri2 = 0. Equations (12) and (13) for such rji have only the solution £ti = 'U2 = which defines 



the global minimum manifold Ai^^{ 

or U2 ^ 
imply r]2 



(0,0,0) 
min 



{x G Ai \ Ui = U2 = 0}. Now let ?73 7^ and ui ^ 
or -112 7^ 0. In this case Eqs. (12)-(15) give |r73p = (1 -|- 77i)(l + ri2) and |r73p = riiri2, which 
1 — r/i and jr/sp = 



?7i(l + rii). Then Eqs. (12) and (15) have the solution 



U2 



-Ml 



V3 



ZUi, 



V2~ 

Vl = -V2 = -Z V2 



(20) 
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where we used the notation z = —(1 + rji)/rj^ G C/{0} and the relation —rj2/Tjl = —z*. 
Note that for a given pair {ui,V2) G C^, z can be any non-zero complex number such 



that {ui, zui, —z*V2,V2) E Ai. The solutions of the form (20) constitute the critical set 
T = {x E Oi\u2 = zui,vi = —z*V2,z G C} C A^saddifi- ^ similar treatment of the region 
O2 shows that the objective function in this region has as critical points only the global 
maximum manifold Ainm/'^ = {x G | = '52 = 0} and the set T. 
Now consider the region O3. 



Let rj^ = 0. Then in the region O3 Eqs. (13) and (14) imply (1 + 772)^2 = rjiVi = 



r]i = 0, i]2 = —1. The solution of Eqs. (12) and (15) for such values of 77, gives the critical 

set {xeM\u^=V2 = 0}cMf^S;Sl 

Let 773 7^ 0. The treatment is similar to the treatment of the case 773 7^ for the region 
Oi and gives the critical set T. A similar treatment of the region O4 shows that the set of 
critical points of the objective function in this region is {x E Ai \ u2 = vi = 0}[JT. 

Combining together the results for the regions Oi, O2, O3, and O4, we find that the 
critical manifolds are the global minimum manifold A4^^^'^\ the global maximum manifold 

A^^ax°\ and the set T[j{x e M\u2 = Vi = 0}[j{x E M\ui = V2 = 0} = M^^^^Sl Since 
4 

[J Oi = Ai, these manifolds are all critical manifolds of the objective function J for the 

2=1 

case w = 0. A simple computation using the constraints (l5])-([7]) shows that the value of the 

objective function at any point x G A^sadd?e equals to 1/2, i.e., J|a^o = 1/2. 

Now we will find Morse indices of the critical manifold Aif^^^l- An arbitrary point 

X = {ui,U2,Vi,V2) G -Msaddk! ^au bc moved into the point x = {ui,U2,Vi,V2) G A<sadd°e ^^^h 
Ul = 0, V2 = by the following transformation: 

Ul = aui + (3u2, U2 = —(3*ui + Q*U2, 
Vl = avi + /9f2, V2 = -/3*vi + a*V2, 

where a,/3 G C, |ap + = 1. For example, a = —^z for x = {ui,zui, —z*V2,V2) G T. 

As in the analysis of the Morse indices for the case 2, in some neighborhood of x we can 
introduce the coordinates ipi, tpi, ^k, yi,p, q, where i,j = 1, . . . , 7 and k,l = 1, . . . , 6. These 
coordinates satisfy the constraint: 



'■^i - If 

i=l / \ i=l 

The second differential of J in these coordinates has the form: 

6 6 



dV = ^dx,2-^dy,^ + 0-|dpr- (21) 



i=l i=l 



It is easy to see that the tangent space to saddle point x is spanned by the vectors 

d d d d 



d^pi dipi dRep dlmp 
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Therefore A^gaddie nondegenerate, dimM.f^^^l = 16 and the Morse indices of A^gaddie 
Z/+ = z/_ = 6. 

Case 3. ||w|| = 1. In this case A_ = 0, A+ = 1, and 

J[Ul,U2,Vl,V2] = ll^tif • (22) 

Let Ui = {x E Ai\vi ^ 0}. Clearly, points in Ui with Ui = form the global minimum of 
the objective. Assume that there are critical points in Ui such that ui ^ 0. For such points 



Eqs. (12)-(15) imply the system of equations 



1^3 1 = V1V2 (23) 
|r/3|' = V2{l + Vi) (24) 

which has only the solutions with t]2 = rj^ = 0. But in the region Ui, vi ^ and therefore 



Eq. (14) implies t]i = 0. Then, Eq. (12) for r]i = r]2 = 1]^ = has the solution -01 = 
which contradicts the assumption ui 7^ 0. As a result, the only critical points in Ui are with 
Ui = 0. These points form the global minimum manifold A^mm = {x E Ai \ui = 0}. This 
manifold is diffeomorphic to the space bundle with 5*^ as a base and S^^ as a fibre. Thus, 
A^^in is connected. We can use ui as normal coordinates in some neighborhood of A^min- 
Thus Aimin is nondegenerate. 

The treatment of the region U2 = {x E Ai\ui 0} is equivalent to the previous 
consideration. The critical points in this region form the global maximum manifold Aimax ~ 
{x E Ai\vi = 0}. Note that U1UU2 = Al. Therefore, all critical points of J correspond to 
the global minimum J = and global maximum J = 1. The critical manifolds corresponding 
to the minimum and the maximum are connected and nondegenerate. □ 

Remark 1 The critical manifolds in terms of the original parametrization of the Kraus 
operators by (mi, U2, "fi, ^2) can be obtained by expressing Ui and Vi in terms of Ui and f j. If 
7^ 0, then it follows from ([8| and (|9|) that 



Ui = jlUi + -j :l'U2, U2 = —VUi + -j -MU2 

, Zq _ Zq 

Vl = fIVi + 7UV2, V2 = -VVi + -JiV2 

\Zq\ \Zq\ 

Thus, for 2:0 7^ and < ||w|| < 1 the critical manifolds are the following: the global 
minimum Aimin = {a; G Ai | Mi = M2 = 0}, the global maximum A^^ax = {x E M.\vi = 
V2 = 0}, and the saddles Ai^ = {x E Ai\u2 = z±Ui, Vi = —z^V2}- Here z± = 2;q/(7 ± ||w||). 
For Zq ^ and ||w|| = 1 (hence 7 7^ 1), the critical manifolds are Aimin = {x E M.\u2 = 
ZqUi/{,1 ~" 1)}; Almax = {x E M.\v2 = ZqVi/{;~^ ~ 1)}) ciud there are no saddles. 

If zo = and 7 > 0, then ui = ui, U2 = U2, Vi = Vi, and V2 = f2- Thus for 7 = 0, 
■^min °^ = {x E M\ui = U2 = 0}, AlL^ax'"^ = {x E M \ Vl = V2 = 0} , and Al^X^ = {^^ 
M\u2 = zui,Vi = —z*V2,z E C}\J{x E M\ui = V2 = 0}. For < 7 < 1 the critical 
manifolds are M^^^''^^ = {x E M \ ui = U2 = 0}, Almax'^'' = {x E M \ vi = V2 = 0}, and the 
saddles AlL°'°'^^ = {x E M\ui = V2 = 0} and Mf^'"'^ = {x E M \ u2 = Vi = 0}. For 7 = 1, 
-^m^n = {xeM\ui = 0} and AlL?/^ = {xeM\vi=0}. 
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If 2;o = and 7 < 0, then ui = U2, U2 = mi, Vi = V2, and V2 = fi- In this case for 
— 1 < 7 < the critical manifolds are the following: M.^^^'^^ = {x G A4\ui = U2 = 0}, 

Min'^"'^ = {x e M\vi = v2 = 0}, A^L°'°'^^ = {x e M\u2 = vi = 0} and A^V"'"'"^ = e 
M\ui = V2 = 0}. For 7 = -1, = {x e M\u2 = 0} and M^^^^^'' = {xeM\v2 = 

0}. 

Remark 2 The values of the objective function at the saddle points satisfy the equality 
J+(w) + J_ (w) = 1. This fact is a consequence of the more general symmetry of the objective 
function, defined by the duality map T : A4 Ai such that T{ui, U2, f 1, f 2) = {vi,V2, Ui, U2) 
as J[x; w]+ J[T(a;); w] = 1 for any x G A^. Thus, if the level set Ty^{a) := {x G | J[x, w] = 
a} for some value a G [0, 1] is known then one immediately gets the level set for the value 
1 — a as rw(l — a) = T(rw(a)). 

5 Connectivity of the level sets 

The level set V^{jj) for an admissible objective value /z G [0, 1] is defined as the set of all 
controls x = {ui,U2,vi,V2) G Ai which produce the same outcome value /i for the objective 
function J[mi, M2, fi, ^2; w], i.e., T^{fi) = {x E Ai| J[x;w] = fi} (we omit the subscript 
w in the sequel). In this section it is shown that each level set for the function J[-;w] is 
connected. This means that any pair of solutions in a level set r(/i) is connected via a 
continuous pathway of solutions entirely passing through r(/i). Practically, connectivity of 
the level sets implies the possibility to experimentally locate more desirable solutions via 
continuous variations of the control parameters while maintaining the same value of the 
objective function. The proof of the connectivity of the level sets for the objective functions 
defined by Q is based on generalized Morse theory, which is presented in the remainder of 
this section. Theorem [2] below formulates the conditions for a generalized Morse function 
to have connected level sets. These conditions are satisfied for the objective function w] 
defined by Q, as stated in the end of this section. Formulation of Theorem [2] includes a 
very general class of functions and can be applied to the investigation of connectivity of the 
level sets for situations beyond the scope of this paper, including landscapes for multilevel 
closed and open systems. 

5.1 Connectivity of level sets of generalized Morse functions 

Let M be a smooth compact manifold of dimension d, and let / be a smooth function 
/ : M ^ M. We suppose that the critical set of /, S" := {x G M|d/(x) = 0} is a disjoint 
union of smooth connected sub-manifolds Cj (i = 1, 2, . . . , n) of dimension di. Let /Xj = /|c-. 

For each point x G Cj there exists an open neighborhood U of x and a coordinate system 
{xi} in U such that 

a n f/ = {x G U\xd,+i = ---=Xn = 0}. (25) 
Consider the following matrix 

X G Q. (26) 

l,m=di+l,...,d 



Ji{x) :-- 



d'fix) 



dxidxr, 
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It is easy to see that if {yi} is another coordinate system in U such that 



QnU^iye U\yd^+i 



Vn = 0}, 



(27) 



and 



Ji{x) : 



d'fjx) 
dyidym 



l,m=di+l,...,d 



X e Ci 



(28) 



then 



rankJj(x) = rankJj(a;). 



(29) 



Therefore we can give the following 

Definition 2 The point x & Ci is said to be nondegenerate z/det Ji{x) ^ 0. 

Definition 3 A critical sub-manifold Cj is said to be nondegenerate if Wx e C^, x is a 
nondegenerate point. 

Let X G Ci and \f{x), \^{x) be the numbers of positive and negative eigenvalues of the 
matrix Ji{x). It is clear that X'l{x), \i{x) do not depend on the choice of coordinate system 
{xi} in the neighborhood of x. One can prove that Xf{x) and X^ix) do not depend on the 
point X & Ci {Xf{x) and X~{x) are continuous and Cj is connected.). Let := X^{x) and 
A~ := Aj~(x). A^ and A^" are called the indices of Cj. 

Definition 4 Let M be a smooth compact connected manifold and / : M — >■ M. Suppose that 
the critical set of f is a disjoint union of (compact) connected nondegenerate sub-manifolds 
Ci. In this case we say that f is a generalized Morse function. Sub-manifolds Ci are called 
the critical sub-manifolds of f. 

Theorem 2 Let M be a smooth compact connected manifold and f be a generalized Morse 
function. Let Ci, i — 1, . . . ,n be critical sub-manifolds of f and iXi — We can assume 

that //min := yWl < < . . . < /Xn =: Aimax- SuppOSC that the sub-manifold Cmax := /"^(Atmax) 

is connected. Suppose also that \/i = l,...,n — 1 the indices > 2, A^ > 2. Then 
V/x : A*min < IJ' < A*max thc sct T^fi) := f~^{ij) is connected. 

Proof. We decompose the proof of the theorem into a sequence of several Lemmas. 

Lemma 1 There exists an open neighborhood U of Cmax such that U is diffeomorphic to 
some bundle E with the base Cmax and the fibre Bd-d„- Here is a k-dimensional ball. 

Proof. M is a compact. Therefore there exists a Riemann metric g G sym(T*M ® T*M). 
(Here T*M is a cotangent bundle of M.) By definition, >C is a restriction of the tangent 
bundle TM to Cmax- Let jV be a sub- bundle of C such that Vx G Cmax the fiber jV^. of J\f 
over X is a subspace of T^M consisting of all vectors orthogonal to T-^Cmax- Let Bi be a 
sub-bundle of Af such that \/x G Cmax the fiber {Bi)x of Bi is a sct of all vectors v of 
satisfying the following inequality: \\v\\ < I (with respect to the metric g). 
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Let ■jy{x){t) {x e M, V e TxM, t e M) be a geodesic line, i.e., the solution of the following 
ordinary differential equation 

V^4.)(t)7.(^)W = (30) 
with the following initial conditions 

7^,(a;)(0) = X, 
%{x){t)\t=o = V. (31) 

Here is a Levi-Civita connection on M with respect to the metric g. The solution of this 
differential equation is defined on the whole real line because M is compact. 

Let Fi for / G (0, +oo) be a map Bi ^ M which assigns to each point {x, v) G Bi 
{x G Cmax,'i^ G {^i)x) the point 7„(a;)(l). It follows from the inverse function theorem that 
there exits a number Zo > such that Fi is a diffeomorphism on its image for alH : < Z < Zq- 
□ 

Lemma 2 If e is small enough then W/i : /^max > IJ> > A*max — £ the set r(/x) = f~^{n) is 
connected. 

Proof. Let Iq be a number from the previous Lemma. It follows from the Morse Lemma that 
for every x G Cmax we can choose coordinates zi, . . . , z^-dn {.^io)x in some neighborhood 
U of zero such that 

foFi,\u^zl + ... + zl_,^. (32) 

Moreover, from construction of these coordinates it follows that in some neighborhood of 
every point xq G Cmax they are differentiable functions of x. Therefore, there exists a 
finite covering {t/i}i=i,...,g of Cmax by open connected sets and a family of diffeomorphisms 
gi : Ui X Bd-dn ~^ '^~^iUi) (i = 1, . . . , g) on its image commuting with the projections such 
that 

/ o ^io ° = + ■ ■ ■ + ^d-dn^ i^l,...,q. (33) 

Here Zi, i — 1, . . . ,d — dn are some coordinates on the ball Bd-d„ and tt is a canonical 
projection from B;g to Cmax- 

We now prove that for every /i : < /i < /q there exists ei > such that V/x : /Xmax — < 
A* ^ A*max, r(A*) C Fi^Bi^. Suppose that Vn = 1,2,... there exists a point a;„ such that 
f{xn) > A*max — l/u and Xn ^ Fi-^^Bi^. Because M \ Fi^Bi^ is compact, then there exists a 
point Xo G M \ Fi^Bii and sub-sequence {xn^.} of {x^} such that x^^, Xq as k ^ oo. We 
find that f{xo) — //max and xq G Cmax- This contradiction proves our statement. If h is 
small enough then Bi-^ Dili C giiUi x Bd-d„) for all i = 1, . . . , g. Therefore if /i > //max — £i 
then f~^{[j) r\Ti~^{Ui) C gi{Ui x Bd~dn) connected. So we find that f~^{lJ>) is connected 

if /X > Aimax-£l- □ 

Lemma 3 Suppose that for some /i : /li < fj, < /Xj+i (i = 1,. . . ,n — 1) the set r(//) is 
connected. Then V// such that /li < /i < /li+i, the set r(//) is connected. 
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Proof. Let u E M. : fii < u < /Xj+i. Let us prove that r(z/) is connected. We can assume 
that u < and let £ be a positive number such that fj,i<u — e<fj, + e< /Xj+i. Consider 
the following sets 

Us = {x\iy — e < f{x) < fi + e}, 

Ue = {x\u-e < fix) < fi + e} (34) 
Consider also the following differential equation on M 

grad/(7(t)) 
= ||grad/(o<0)ir 

(Recall that M has a Riemann metric). The right hand side of this equation is well defined 



on Us- The solution of (35)is 7x(^) with the initial condition 

7a;(0) =x, X e T{fi). (36) 

By the extension theorem [l9l [50] this solution must leave the compact set f/e/2- It is easy 
to prove that fijxif)) = t + fi. So the solution 'jxit) is defined and unique on the interval 
iy — H — £^/3, /i + e/3). Therefore we have a smooth map /S.^^y : T{n) r(z/), x ^— --yxiy — /w)- 
By the same means we can construct the map A,^^^ : T{u) r(/i). A^j^(a;) = y ii and only if 



a; and y lie on the same integral curve of (35). We have A^^^^o A,^^^ = id and A,^^^o A^_,^ = id. 
So r(/x) and r(i/) are diffeomorphic. □ 

Lemma 4 Suppose that the assumptions of the theorem hold. Let /i G M.' /ij < /i < yUi+i, 
fii = 2, . . . ,n — 1, and r(/i) zs connected. Then Vz/ swc/i i/iai yUj_i < u < fii, the set r(z/) zs 
a/so connected. 



Proof. We prove this lemma only for the case of connected Cj. The general case is analogous 
to this case. 

As in Lemma [T| let Bi be a bundle with the base Cj which consists of all vectors v 
normal to Ci and such that ||f || < 1. We have Bi-^ C Bi^ for li < I2. Let Fi be a map 
Bi ^ M constructed as in Lemma 1. As in Lemma [ij we find that Fi is a diffeomorphism if 
< / < /o for some positive number /q. As in Lemma 1 we find that for every I'q < Iq there 
exists a covering {t/j}j=i,...,p of Ci by open connected sets and the family of diffeomorphisms 
gj : Uj X Bd-di — > n~^{Uj) on its image commuting with the projections such that 

foFiogj = zl + ... + zl+- zl+ - ... - z^d-d, + /ii (37) 

Here Bd-di is a d — rfj- dimensional ball and tt a canonical projection from Bi' to Cj. 

It is easy to see that for every I'q < Iq there exists a positive number li < I'q such that 
Vj = 1, . . . ,p Bi^ f\Ti^^{JJj) C gjiJJj X Bd-di)- every /i < I'q there exists a positive number 
62 such that Vx G Cj, (i3;^/2)a; n F[^^{V{iJii + /t)) 7^ Vk : |k| < £2- We now prove that 
Bi^ U TX^^(Uj) n F^~"'^(r(/i + k)) is connected Vj = 1, . . . ,p if |k| < 62- Indeed, let xi and X2 
be two points which lie in the set Bi^ U Ti~^{Uj) fl F^~^(r(yU + n)). We can consider only the 
case K > 0. The set gj{Uj x 5d_dJ n F^~^(r(/ij + k)) is diffeomorphic to R-^^ x S'^r-i x Ui 
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and connected. Let 'j{t) t G [0, 1] be a path in gj{Uj x Bd-di) H F^^^(r(/ij + k)) such that 
7(0) = Xi, 7(1) = X2- Let d{x) be a function on defined as follows: d{{z,v)) = ||f||^, 
where z & Ci and G {Biq)^- Let x G Fj^"'^(r(/i + k)) fl x Bd^dJ and be a 

projection of Vd{x) to the tangent space of F;~^(r(/i + k)) at x. It is obvious that w{x) 7^ 
Vx G Fj7^(r(/i + k)) n gj{Uj X Bd-di) ^ i^i'o \ ^h) if ^0 ^ sufficiently small number. So 
we can retract the path 7(t) along the vector field w to the part 7(t) which lies in Bi-^ and 
connects the points xi and X2. So U 7r~-'^(f/j) fl Fj~^(r(/i + k)) is connected. Now we can 
find that Bi-^ U nF;~"^(r(/i + /«)) is connected. 

Now let xi,X2 G r(/i), /i < /ij, \fi - < 63. Let U = Fi^{B)i^/2, V = Fi^{B)ij3, 
W = Fi^^{B) 1-^/4. At first suppose that xi ^ U and X2 ^ U. Let 7x1 (^), 7^2 (^) be solutions of 



the differential equation (35) with initial conditions Xi and X2 respectively. The paths 7^1 (^) 
and 7x2 (^) intersect the sub-manifold T^fi + es) at the points yi and 1/2 if £3 is enough small. 
Let 6{t),t G [0, 1] be a path such that Vt G [0, 1] 6{t) G r(/i + e) and yi = 5(0), ?/2 = 5(1). 
We must consider the following two cases. 

1) 5 n y = 0. If £3 is small enough then we can deform the part 6 along the vector field 
V//||V/|p to the part 6 which lies on T{^) and connects the points xi and X2- 

2) 6 n V ^ ^. If £3 is small enough then yi,y2 ^ V. We can decompose the part S as 
6 = ai o p o 0(2, where 

5^(1) g9\/, yt e[0,l]^{t) 

51(0) G dV, yt G [0, 1] SI(t) ^ V. (38) 

If £3 is a sufficiently small positive number we can deform the paths 5i and 02 along the 
vector field V//||V/|p into the the paths ai,a2 C r(/i) such that ai,a2 ^ and a2(0) = 
Xi, 0:2(1) G t/, ai(l) = X2 and ai(0) G f/. But, it has been proved that UnT{fi) is connected. 
Therefore, there exists a path /? C r(/i) such that /5(1) = q;i(0) and 02(1) = /?(0). We see 
that the path ai o f3 o a2 connects the point Xi and X2- 

Consideration of the case with xi E U or X2 E Bu is analogous to consideration of the 
previous case. The statement of the theorem follows from these four Lemmas. □ 

Theorem 3 Each level set of the objective function J[-, w] defined by ^ is connected. 

Proof. The objective function J[-,w] is a generalized Morse function. The sub-manifold 
of solutions corresponding to the global maximum in the coordinates Ui,U2 ^ is defined 
by IImiII = ||'U2|| = 1, {ui,U2) = 0. It is a Stiefel manifold, A^^ax = V2(C^), and hence 
is connected. The Morse indices of the function J[-,w] are zz-t > 2 at any saddle sub- 
manifold. Therefore this function satisfies the conditions of Theorem |2] and its each level set 
is connected. □ 



6 Conclusions 

In this paper the landscape of the objective functions for open quantum systems controlled 
by general Kraus maps is investigated in detail for the two-level case. It is shown that 
a typical objective function has: (a) no false traps, (b) multi-dimensional sub-manifolds 
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of the optimal global solutions, and (c) each level set is connected. These results may 
be generalized to systems of arbitrary dimension A^, although a full enumeration of the 
critical sub-manifold dimensions remains open for analysis. The landscape analysis and the 
conclusions rest on assuming that the controls can manage the system and the environment. 
Managing the environment, in practice, is likely not highly demanding, as control over only 
the immediate environment of the system is most likely needed. The critical point topology 
of general controlled open system dynamics could provide a basis to explain the relative 
ease of practical searches for optimal solutions in the laboratory, even in the presence of an 
environment. 
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